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Abstract. The classical dynamical model for reactions induced by weakly-bound 
' nuclei at near-barrier energies is developed further. It allows a quantitative study 

of the role and importance of incomplete fusion dynamics in asymptotic observables, 
such as the population of high-spin states in reaction products as well as the angular 
distribution of direct alpha-production. Model calculations indicate that incomplete 
^ , fusion is an effective mechanism for populating high-spin states, and its contribution 

(— I ' to the direct alpha production yield diminishes with decreasing energy towards the 

Coulomb barrier. It also becomes notably separated in angles from the contribution 
of no-capture breakup events. This should facilitate the experimental disentanglement 
^ ' of these competing reaction processes. 

in : PACS numbers: 25.60.Pj, 25.60.Gc, 25.60.-t, 24.10.-i 
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Q I 1. Introduction 

Nuclear physics research has entered a new era with developments of radioactive nuclear 
^ . beam facilities, where nuclear reactions are the primary probe of the new physics, such 
■ as novel structural changes. In those facilities, the low-energy nuclear reactions research 
is highly focused on understanding astrophysically important reaction rates involving 
exotic nuclei. These are often weakly-bound with a few-body, cluster structure that 
can easily be dissociated in their interaction with other nuclei. Understanding the 
breakup mechanism and its impact on nuclear reaction dynamics is essential. A major 
consequence of breakup is that a rich scenario of reaction pathways arises, such as 
events where (i) not all the resulting breakup fragments might be captured by the target, 
termed incomplete fusion (icf), (ii) the entire projectile is captured by the target, called 
complete fusion (cf), and (iii) none of the breakup fragments are captured, termed no- 
capture breakup (ncbu). 

Since the availability of intense exotic beams is still limited, extensive experimental 
research has recently been carried out exploiting intense beams of stable weakly-bound 
nuclei, such as ^'^Li and ^Be [H E]. Understanding the effect of their breakup on near- 
barrier fusion has been a key aspect of these investigations |3]. These have definitively 
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demonstrated that breakup suppresses the above-barrier fusion cross sections. Most 
recently, experimental activities have been focused on disentangling breakup and 
competing reaction mechanisms from inclusive and exclusive coincidence measurements 
[H El El [TJ [8] . One of the challenges is to obtain a complete quantitative understanding 
of the breakup mechanism and its relationship with near-barrier fusion. This research 
is guided by complete sub-barrier breakup measurements p]. 

Theoretical works have addressed the low-energy reaction dynamics of weakly- 
bound nuclei using quantum mechanical, classical and mixed quantum-classical 
approaches [TOl [TTl [121 [131 [IH [IHl [HI [17]. Among these, the continuum- discretised 
coupled channels (CDCC) framework has been very successful [ini [IHl [lEl [IT] . However, 
existing quantum models have limitations [18], as they cannot calculate integrated ICF 
and CF cross sections unambiguously. Neither, after the formation of ICF products, 
can they follow the evolution of the surviving breakup fragment (s) since ICF results in 
depletion of the total few-body wave function. 

These difficulties are overcome by the three-dimensional classical dynamical 
reaction model suggested in Ref. [19]. A crucial input of this model is a stochastically 
sampled breakup function proposed in Ref. [20], which can be determined from sub- 
barrier breakup measurements [9l [20] . This function encodes the effects of the Coulomb 
and nuclear interactions that cause the projectile breakup. Hence, this approach is not 
a breakup model, rather it is a quantitative dynamical model for relating the sub-barrier 
NCBU to the above-barrier ICF and CF of weakly-bound nuclei [20] . The model has been 
successfully applied to interpreting fusion measurements of weakly-bound nuclei [3 [9] , 
and isomer ratio measurements [21j . 

In this paper, I report on further developments of this classical dynamical model. 
The key new aspect is the time propagation of the surviving breakup fragment and 
the ICF product, allowing the description of their asymptotic angular distribution 
and dynamical variables. These should be very useful in (i) current experimental 
activities aimed at disentangling competing reaction mechanisms from asymptotic 
observables such as alpha-production yields [8], (ii) interpreting particle-7-coincidence 
experiments [22], and (iii) applications to 7 ray spectroscopy [231 1211 [25]. The new 
developments are illustrated with a simplified test case that does not aim at adjusting 
any measurements. It is worth mentioning that various dynamical models have been 
proposed for multi-fragmentation and ICF in heavy-ion induced reactions at energies 
well-above the Coulomb barrier (> 10 MeV/nucleon), as recently reviewed in Ref. [22] . 
In this different context, some of those models (e.g., Refs. [261 l2Z]) use concepts and 
techniques (e.g., classical trajectory, Monte Carlo sampling) that overlap with some 
involved in the present model. The model is explained in Section 2. In Section 3, 
numerical results are discussed, and a summary is given in Section 4. 

2. Model 

The main features of the model are as follows: 
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(i) The target T is initially at rest in the origin of the laboratory frame, whilst the 
weakly-bound (two-body) projectile P approaches the target (along the z-axis) 
with incident energy Eq and orbital angular momentum Lq. For each Lq (chosen 
to be an integer number of h) an ensemble of incident projectiles is considered. 
Including the P — T mutual Coulomb and nuclear forces, classical equations of 
motion determine an orbit with a definite distance of closest approach Rmin{Eo, Lq). 

(ii) The complexity of the projectile dissociation is empirically encoded in a density of 
(local breakup) probability V^u{R), a function of the projectile-target separation 
R, such that V^jj{R)dR is the probability of breakup in the interval R to R-\- dR 
(see Appendix A). A key feature is that for a given projectile-target combination, 
both measurements [9l [20] and CDCC calculations [19] indicate that the integral 
of this breakup probability density along a given classical orbit is an exponential 
function of its distance of closest approach, Rmin{EQ, Lq): 



Consequently, V^ui.R) has the same exponential form, Vgu^R) oc exp(— ai?). 
(The factor of 2 highlights that breakup may occur along the entrance or exit 
branch of the trajectory, although the exponential form will clearly place maximum 
probability of breakup at Rmin-) This function is sampled to determine the 
position of breakup in the orbit discussed in (i). In this position, the projectile 
is instantaneously broken up into fragments Fl and F2. These interact with T, and 
with each other, through real central two-body potentials having Coulomb barriers 
y^"' at separations i?^, i,j = 1,2,T, i ^ j. 

(iii) The instantaneous dynamical variables of the excited projectile at breakup, namely 
its total internal energy eu, its angular momentum £12 and the separation of the 
fragments di2 are all Monte Carlo sampled. The initial separation di2 between 
the fragments and its orientation du are determined by sampling the radial and 
angular probability distributions of the projectile ground-state (g.s.) wave function, 
respectively. For a two-body projectile with 0+ g.s., a very good approximation for 
calculating du is through a Gaussian sampling function in the classically allowed 
region of the fragments, whilst the orientation of di2 is isotropic. (This will be used 
in the test case below.) For high £12 excitations, when there is no barrier between 
Fl and F2, di2 is equated with their external turning point. The orientation of £12 is 
chosen randomly from all directions orthogonal to d^. iu is sampled uniformly in 
the interval [0,£m.aa]; whilst for £12 an exponentially decreasing function for energies 
between the top of the barrier (V^^) and a chosen maximum Emax is sampled. Both 
imax and Emax are increased until convergence of the observables occur. 

(iv) Having fixed the position and dynamical variables of the excited projectile 
fragments at the moment of breakup, the instantaneous velocity of the particles Fl, 
F2 and T is determined by conservation of energy, linear momentum and angular 
momentum in the overall center-of-mass frame (see Appendix B). These breakup 




(1) 
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initial conditions are transformed to the laboratory frame where the three bodies 
are propagated in time. The calculated trajectories of Fl, F2 and T determine the 
number of ICF, CF and NCBU events, fragment Fj being assumed to be captured if 
the classical trajectories take it within the fragment-target barrier radius -R]J. 

(v) From the N breakup events sampled for each projectile angular momentum Lq, 
the numbers of events Ni in which i = (ncbu), 1 (iCP), or 2 (cp) fragments are 
captured determine the relative yields Pi = Ni/N of these three reaction processes 
after breakup, with Pq + -Pi + A = 1- The absolute probabilities Pi^Eq, Lq) of these 
processes are expressed in terms of the relative yields and the integrated breakup 
probability over the whole trajectory PBu{Rmin)- 

Poi^O, Lq) = PBuiRmin) Pq, (2) 

PlipQ, Lq) = PBu{Rmin) Pi, (3) 
P2{Eq, Lq) = [1 — PsuiRmin)] H{Lcr — -^^o) 

+ PBuiRmin) P2, (4) 

where H{x) is the Heaviside step function and Lcr is the critical partial wave for 
projectile fusion. The cross sections are calculated using 

(T,(Eo) = ttA^ ^(2Lo + l)P.(^o, ^0), (5) 

Lo 

where = /i^/[2mpi?o] and mp is the projectile mass. 

Beside the absolute cross sections ([5]), asymptotic observables, such as the angle, 
kinetic energy and relative energy distributions of the fragments from ncbu events, 
are calculated by tracking their trajectories to a large distance from the target. 

(vi) For the ICF events, the time propagation of the ICF product and the surviving 
breakup fragment is now incorporated into this model as follows. The fragment Fj 
(after overcoming the Coulomb barrier V^^) reaches the target radius, forming the 
ICF product, while the other fragment flies away. At this moment, the three-body 
propagation turns into a two-body propagation, with definite interaction potentials 
and initial conditions. These are determined by the position and velocity of the 
three particles, at the moment when the ICF product is formed. These also yield 
the spin and excitation energy distributions of the primary ICF product. The 
asymptotic angular distribution of the ICF product and the surviving breakup 
fragment is calculated in terms of their trajectories. 

3. Numerical results 

This model is implemented in the platypus code, as described in Ref. [28]. In order 
to illustrate the new developments of (vi), calculations are carried out at a laboratory 
energy of Eq = 45 and 65 MeV, for the reaction of a pseudo-®Be projectile P (modeled 
as a weakly-bound s-state of two a-particles [T9|) with the ^°^Pb target T. 

The breakup function PBu{Rmin) in eq. ([1]) has parameters A = 5.98 x 10^ and 
a = 0.85 fm~^, which was deduced in Ref. [19J from mapping fusion measurements for 
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Figure 1. Angular distribution of direct alpha-production for *Be -I- ^"^Pb at two 
laboratory energies Eq: (a) 65 MeV, and (b) 45 MeV. With decreasing energy towards 
the Coulomb barrier, the NCBU events dominate, separating its centroid substantially 
from that of the ICF events. The total alpha-production distribution changes its shape 
notably. 

the ^Be -|- ^ospi^ system with this test reaction. This is because the model is hmited to 
solving a three-body problem. However, very recent measurements [9] have shown that 
prompt ^Be breakup occurs dominantly through an excited ^Be nucleus, validating the 
approximation of a ^Be projectile by ^Be. 

The nuclear interaction between the alpha particle and the ICF product ^^^Po is 
considered to be the Woods-Saxon potential {V, r, a) = (33.98 MeV, 1.48 fm, 0.63 fm) 
resulted from the global Broglia-Winther parametrization [29] . (Please note that in the 

1 /3 

potential the radius parameter is multiplied by .) The rest of the model parameters 
are the same as in Ref. [19j . 

Figure [1] shows the angular distribution of direct alpha-production for two 
laboratory energies near the P — T s-wave Coulomb barrier (39.9 MeV), namely (a) 
Eq = 65, and (b) 45 MeV. The contribution of the ICF and NCBU events is represented 
by thick solid and thick dashed lines, respectively. Their incoherent sum is represented 
by the thin solid line. Its shape significantly changes as the incident energy decreases. 
While the contribution of the ICF and NCBU events appears to be similar at well-above 
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Figure 2. ICF and ncbu cross sections as a function of the relative angular momenta 
Lq for ^Be + ^°*Pb at Eq = 45 MeV. The contribution of high-partial waves shifts the 
NCBU distributions in Fig. [T]to smaller angles, with respect to the ICF distributions. 



barrier energies [panel (a)], the ncbu contribution gradually dominates with decreasing 
energy towards the barrier [panel (b)]. Here, its centroid significantly separates from the 
centroid of the ICF contribution. Both centroids shift to higher angles as the incident 
energy decreases, due to the reduction of relative partial waves affecting these reaction 
processes. However, the ncbu centroid always remains lower than the ICF centroid, as 
expected. This is because higher partial waves contribute to the ncbu process (see Fig. 
E]). Fig. [2] presents the incident angular momentum distribution of the ICF (solid line) 
and NCBU (dashed line) processes at a laboratory energy of = 45 MeV. 

Figure |3] shows the angular momentum distribution of primary ICF (^^^Po) and 
CF (^^'^Rn) products at (a) Eq = 65, and (b) 45 MeV. The ICF spin distribution is 
represented in terms of both the ^Be incident angular momentum Lq (thick solid line) 
and the angular momentum brought in by the a-particle into the ^osp]^ target (thin 
solid line). These representations result in two very different ICF angular momentum 
distributions, both providing relevant features of the ICF dynamics as explained below. 

Comparing the thick solid line with the CF spin distribution (dashed line) it is 
seen that the localization of the Lo-window for ICF significantly depends on the incident 
energy Eq. Their overlap strongly increases as Eo decreases towards the barrier [panel 
(b)], indicating that ICF and CF are two competing reaction processes at near-barrier 
energies. However, CF following ^Be breakup is here a very small component (1.5%) of 
the CF cross section (149.2 mb). Instead, either one of the a-particles is captured by the 
target, contributing to the ICF cross section (84.7 mb), or the two a-particles survive and 
contribute to the ncbu cross section (238.8 mb). At Eq = 65 [panel (a)], CF following 
'^Be breakup is very substantial, representing 41% of the CF cross section (981.5 mb). 
This significantly shifts down the maximum of the CF spin distribution (dashed line) 
with respect to the critical angular momentum L^r for ^Be fusion as an inert projectile 
(thick arrow). Here, Lq values around L^r determine the ICF cross section (273.6 mb) 
which is similar to the ncbu cross section (259.3 mb). 
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Figure 3. Spin distribution of primary fusion products resulting from the ^Be + 
208pj-, reaction at two laboratory energies Eq: (a) 65 MeV, and (b) 45 MeV. The 
CF (^^^Rn) spin distribution is presented by the dashed line, whilst the ICF (^^^Po) 
spin distribution is shown in two representations, i.e., in terms of (i) the '^Be incident 
angular momentum Lq (thick solid line), and (ii) the angular momentum brought in 
by the a-particle into the ^''^Pb target (thin solid line). The thick arrow denotes 
Lcr for ^Be fusion as an inert projectile, whereas the thin arrow is Icr for fusion of 
a direct beam of a-particles (with half of Eq) on ^"^Pb. The relevant features are: 
(1) the ICF mechanism increases the angular momentum brought in by the a-particle 
into ^°*Pb (thin solid line goes beyond the thin arrow), (2) the relative localization 
of the Lg-window for CF and ICF significantly depends on Eq (dashed and thick solid 
lines), strongly overlapping as Eq decreases, and (3) CF following ^Be breakup becomes 
substantial with increasing Eq, shifting down the maximum of the CF spin distribution 
(dashed line) with respect to the thick arrow. 

The ICF spin distributions (thin sohd hne) in Fig. |3] are the crucial ones in order to 
access the effectiveness of the ICF mechanism for populating high-spin states in ^^^Po. 
It is observed that the tail of these distributions goes well beyond the thin arrow which 
denotes lcr for fusion, on the ^osp]^ target, of a direct beam of a-particles with half of 
Eq. Clearly, the extra torque caused by the interaction between the fusing and surviving 
a-particles enhances the angular momentum of the ICF product. This interaction also 
affects the excitation energy of ^^^Po substantially, as presented in Fig. HJ 

Fig. m shows the excitation energy distribution of the primary ICF product ^^^Po 
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Figure 4. Excitation energy distribution of the primary ICF product ^^^Po resulting 
from the ^Be + ^°*Pb reaction at two laboratory energies Eq: (a) 65 MeV, and (b) 45 
MeV. The arrow denotes the excitation energy of ^^^Po, when this compound nucleus 
is formed by fusion of a direct beam of a-particles, with half of Eo, on ^°*Pb. Following 
breakup of ^Be, the fusing a-particle is significantly affected by the interaction with 
the surviving a-particle. This spreads the ^^^Po excitation energy over a range of 
values. The distribution becomes narrower and its maximum approaches the arrow, as 
Eq decreases towards the barrier. 

(thin solid line) at (a) Eq = 65, and (b) 45 MeV. It is caused by the interaction between 
the fusing and surviving a-particles during the ICF process. The arrow corresponds to 
the pre-determined excitation energy of ^^^Po, when this nucleus is formed through CF of 
a direct beam of a-particles (with half of Eq) on ^°^Pb. This value remains smaller (by a 
few MeV) than the excitation energy associated with the maximum of the distribution. 
The shape of the distribution also changes with Eq, becoming narrower as Eq decreases 
towards the barrier. These distributions along with the ICF spin distributions in Fig. [3] 
(thin solid line) are vital for a reliable prediction of the yield of final ICF products. To 
my knowledge, these aspects have not been yet included in broadly used evaporation 
codes which mostly describe CF evaporation residues. 

The present model can be applied to more realistic cases (e.g., ^^Be induced 
reactions), provided all the necessary information for Monte Carlo sampling is known 
(e.g., the breakup function and the projectile g.s. wave-function). When one of the 
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breakup fragments is neutral, its capture by the target nucleus can be assumed to occur 
when its trajectory takes it within the target radius. Although not implemented yet, the 
model can also predict the elastic angular distribution of the projectile. This distribution 
can easily be calculated in terms of the projectile-target orbits and the probability of the 
weakly-bound projectile's survival (see Appendix A), i.e., this probability weighting the 
contribution of the different projectile-target orbits. Direct reaction processes other than 
the elastic breakup may significantly contribute to observables associated with NCBU, CF 
and ICF. These are not included in the model yet. However, their contributions could be 
included using transfer and inelastic-breakup functions. Of course, these would make 
the calculations more complex, as these direct reaction processes result in additional 
bifurcation points along the projectile-target orbits. Nevertheless, their inclusion as 
well as the treatment of more complex projectiles (e.g., ®He and ^^Li) will be interesting, 
necessary developments of the present approach. 

4. Summary 

The classical dynamical model for reactions of weakly-bound nuclei at near-barrier 
energies has been developed further. It allows us to quantify the role and importance 
of ICF dynamics in asymptotic reaction observables, such as the angular distribution 
of the direct alpha-production. The ICF contribution to this yield diminishes with 
decreasing energy towards the barrier. However, the ICF and NCBU contributions 
are clearly separated in angles at near-barrier energies, facilitating the experimental 
disentanglement of these competing reaction processes. The present developments also 
provide the spin and excitation energy distributions of primary ICF products, which 
are essential ingredients for calculating the yield of final ICF products with present 
evaporation codes. The ICF mechanism appears to be an effective route for producing 
high-spin states. All these observables may also be affected by other direct processes, 
such as transfer [30], which are not included in the model yet. Nevertheless, the 
present classical dynamical model is a powerful tool for interpreting fusion measurements 
involving radioactive nuclei and in applications to 7 ray spectroscopy. The development 
of a unified quantum dynamical description of relevant reaction processes (iCF, CF, NCBU 
and transfer) remains a great theoretical challenge. One possibility of tackling this issue 
could be through a time-dependent density-matrix approach incorporating the concept 
of quantum decoherence [HI [31] . 

Acknowledgments 

The author thanks B. Dayman and M. Dasgupta for suggesting Appendix A, and J. A. 
Tostevin for discussions related to Appendix D. Support from the UK Science and 
Technology Facilities Council (STFC) Grant No. ST/F012012/1 is acknowledged. 



Modelling incomplete fusion dynamics of weakly-bound nuclei at near-barrier energieslO 



Appendix A. Breakup probability function 

Let us define two probabilities: (i) the probability of breakup between R and R + dR, 
p{R)dR [being p{R) a density of probability], and (ii) the probability of the weakly- 
bound projectile's survival from oo to R, S{R). The survival probability at i? + dR, 
S{R + dR), can be written as follows 

S{R + dR) = S{R) [1 - p{R)dR\. (A.l) 

Expression (1A.1I) suggests the following differential equation for the survival 
probability S{R), 

^ = -S{R)piR), (A.2) 

whose solution is [5'(oo) = 1]: 

S{R) = exp(- f p{R)dR). (A.3) 

J oo 

From (TO]) , the breakup probability at R, B{R) = 1 - S{R). If J^p{R)dR < 1, 
B{R) can be written as 

B(R) ^ I p{R)dR. (A.4) 

J oo 

From (lA.4p . identifying p{R) with Vj^u{R), we obtain expression ([T]) for the breakup 
probability integrated along a given classical orbit. 



Appendix B. Matching prior- and post-breakup stages 

The integrals of motion in the overall center-of-mass (cm) system are the total energy 
~ (mJ+mp) -^0' total linear momentum Ptot = 0, and the total angular 
momentum Lfot = mpbQ{v — Vcm) that is orthogonal to the initial reaction plane. 
tut, mp, bo, V, and Vcm are the mass of the target and projectile, the impact parameter 
between the projectile and the target, the velocity of the incident projectile in the 
laboratory system and the CM velocity, respectively. 

Just after breakup, the two-body projectile is excited to a definite state (ei2, 
and di2), as explained in Section 2. The relative vector between P and T (Rpr) is also 
known. Thus, the separation between the three bodies is known. The modulus of the 
velocity between P and T (Vpt = Ppt/ppt) results from the total energy conservation 

Etot = ei2 + f/ir(rir) + U2T{r2T) + Ppt/^Ppt, (B.l) 

where U is the interaction potential between the target and the breakup fragments. 

The total linear momentum Ptot = Vt + Vi + V2 = Pt + Pp*, where pp* is the 
momentum of the center of mass of excited P relative to the overall CM. We need the 
velocities of P and T relative to the overall CM {vp and vt) to complete the initial 
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conditions for subsequent propagation in time of the three bodies. These velocities are 
related to each other by the expressions 

fT = -—fp, (B.2) 
my 

— * — * — * 

VpT = Vp — Vt, (B-3) 

where the magnitude of VpT is known through expression f IB.ip . To know the direction 
of this velocity the conservation of total angular momentum is applied. 

The total angular momentum Ltot = ^12 + LpT, so the angular momentum {Lpt) 
associated with relative motion of P and T about CM is known. This vector can be 
written as 

— * — * — * 

LpT = mp RpT X Vp. (B.4) 

We now write vp in terms of radial and transverse components as follows: 

Vp = Vp^ f + v^p g, (B.5) 

where r = Rpt/Rpt and q = n x f, being n = Lpt/Lpt- The transverse component 
Vp^ = Lpt/ (mpRpT), and for the target v!^^ = —Lpt/ {niTRpT)- The radial component 
is obtained using expressions flB.2l) - flB.3p and knowing the transverse component: 

»? = ±{v^,-fel?'(l + !^)lT'Al + ^). (B.6) 
I L niT -I mT 

Both positive and negative roots are consistent with the conservation of the integrals 

of motion. Hence, both roots are uniformly sampled. Finally, the position and velocity 

vectors of the projectile fragments and the target are transformed to the laboratory 

system using Galilean transformations. 
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